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Abstract 

We will offer the formal arguments that the very high multiplicity (VHM) events are 
defined by the heavy QCD jets and, with exponential over multiplicity accuracy, the number 
of jets must decrease with the increasing multiplicity. This must rise the mean transverse 
momentum of secondaries with multiplicity. The effectiveness of the LLA in the VHM domain 
is also discussed. 
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CN ■ 1. The interest of experimentalists in the very high multiplicity (VHM) physics has 

O . grown for the last few years But in spite of all efforts there has not been any good 
quantitative theory of such events until now. The reason lies in the special kinematics of 
VHM events: the produced particles have approximately the same small momenta and for 
this reason the usage of LLA becomes problematic. 

The purpose of the present paper is to give the formal arguments that (i) the hadron 
processes become hard in VHM domain and (ii) this tendency to the events hardness is 
! conserved in the asymptotics over multiplicity. 
^1 2. It is useful to introduce the generating function 

o 

p i' where Umax = y/s/m is the maximal multiplicity of hadrons of mass m at given CM energy 

p • \/s. It appears because of the energy-momentum conservation law. We will assume that 

1^ . the total energy y/s can be arbitrary large. The lepton and photon production is neglected. 

' It seems natural to omit the dependence from rimax if n is far from it: in the opposite case 
the dependence on the phase-space boundaries would be sizeable for all n. We will assume 

^ I that z is so small that 

c3 ■ n<^ Umax (2) 

are important in (|l]). Assuming (H), we can consider T{z, s) as an entire function of z in the 
future. 

One may turn ([l|) and define cr„(s) through T{z,s) by the inverse Mellin transform: 
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This integral can be calculated for large values of n by the stationary phase method. The 
corresponding saddle point is defined by the equation: 

d 

n = z—lriTiz^s). (3) 
oz 

Then the asymptotic estimation exists: 

(T„(s) oc e-"^"^=("'^), (4) 
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where Zc{n,s) is the solution of @. The estimation iQj exists if Zc{n,s) > 1, i.e. if n is 
larger than the mean multiplicity n{s); at the same time n can not be too large, see (^. The 
estimation (H) signifies that the quantity 

/i'(n, s) = In an{s)/(ytot{s) (5) 

is defined with the 0(l/n) accuracy only by the solution Zc{n,s). 

As follows from (||), with the increasing multiplicity Zc{n, s) must tend to the singularities 

Zs of T{z, s): 

Zc{n, s) ^ Zs . (6) 

We will assume that |0| T{z, s) is regular in the complex z plane inside the unit circle 
(otherwise cr„ would be the increasing function of n inside the domain (||)). Keeping in mind 
(H), the estimation (H) means the following supposition: 

(I). The asymptotics over multiplicity is governed with the exponential accuracy by the 
leftmost singularity over z. 

It was shown |^ that the multiperipheral kinematics leads to T{z, s) which is regular 
for the arbitrary finite z. For this reason the soft channel of hadrons production leads, as 
follows from (H), to fi'{n,s) ~ In Zc{n,s) increasing with n. On the other hand the pQCD 
jets generating function gives the singularity at Zs{s) — 1 ~ l/nj{s) > 0, where the mean 
multiplicity in the jet is nj{s), Innj(s) oc Vln s 0|. Therefore, for hard processes, 

yu'(n, s) = l/nj{s) + 0(l/n), 

i.e. /i'(n, s) is the n independent quantity which decreases with the increasing energy. 

We conclude that the VHM processes are mostly the hard processes and at first sight 
they can be described by the pQCD. Whether this is right or not is the question of our 
interest. 

3. We will consider particles production in the "deep inelastic scattering" (DIS) kine- 
matics. Let J-'ab{x,q'^;u!) be the generating functional: 



J'abix,q^;uj) = J2 I dn,{k)\{u'^{k1) \aT-''^ {k,, k,, k 

V i=l 
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where a^b^ the production amplitude of u partons (r^ = {q,q,g)) with momenta 

{ki, k2, ku) in the process of scattering of the parton a on the parton b; dVL^ik) is the 
phase space element; uj'^ik'^) is the "probe function", i.e. the correlation functions 

NT-^-^{klkl...X;x,q') = \{jJjg-\ri:Fab{x.q'-,u:) 

The generating functional is normalized on the DIS structure function Vah{x,q^), 

Tab{x, q^; u = 1) = Vab{x, q^), 

which is described in the leading logarithm approximation (LLA) by the ladder diagrams. 
We will consider the approximation when the cutting line passes only through the "steps" 
of the ladder diagram. In this case T>ab{x, q^) has a meaning of the probability to find the 
parton a in the parton b. 
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Considering the VHM region, one must assume that the mass \ki\ of the produced parton 
is large. For instance fcf ^ A^, where A is the virtuality of the parton of the pre-confinement 
phase. Then 

The LLA means that 

A^ « \ql^\ « « \q^\, (8) 

Therefore, \ki\ are also strongly ordered. The approximation (0) means that the produced 
partons longitudinal momenta are small: 

ln{l/xi)<\n\ql,\. (9) 

It is useful to consider the Laplace image over ln(l/a;): 

Tab{x,q^;u;)= [ ^ ( -] F,,(j,g2.^). (10) 



27ri \x ^ 

Then, taking into account the above mentioned conditions, one may find the DGLAP evo- 
lution equation: 

t^Fat{j,t;uj) = J]<^L(j>'WFc6(j,t;^), (11) 

c,r 

where t = Indg^j/A^), 



and -Pac(^) the regular kernels of the Bete-Solpiter equation for pQCD [|]. The equation 
([TT| ) coincides at u;'^ = 1 with the habitual equation for Laplace transform of the struc- 
ture function Vab{x,q^). While the equation ([Tl| ) was being derived, only one additional 
assumption had been used for our problem cu^ = u'^ik'^). 

The dominance of the gluon contributions for the case x -C 1 will be used and for this 
reason all parton indexes will be omitted. One may find the solution of ([TlD in terms of the 
z/-gluon correlation functions N^^\ Omitting the t dependence in the renormalized constant 
Os, let us write: 

F(j, t;^) = D(j,t)exp|^l j \[dU{uo{U) - l)N^^\UM. ■■■^u\x,t)^ . 

In the VHM domain, where x ^ 1 is important, one must consider (j — 1)^1. Then 

Ar«(ti;j,t) = <^(j)~^»l. 

3 - 1 



The second correlator 

N^^\h,t2;j,t) = 

is negligible at (j — 1) ^ 1. Therefore, in the LLA, 

F(j, t; u) = D(j, t) exp |(^(j) ^ dt,{u;{ti) - 1 



Taking u{t) = const, one may find that F{j,t]uj) has the Poisson distribution with the 
"mean multiplicity" ~ ^{j)t- 

If the time necessary to confine a parton into the hadron is ~ {^/^), then the parton of 
the mass (virtuality) ^ A must decay on the partons of lower mass since its life time is 
~ << (1/A). This is the reason of the pQCD jets formation. 

If the quantity 

u{t, z), uj{t, 1) = 1, t = In A;VA^ 
is the generating function of the pre-confinement partons multiplicity distribution 

Qn 

UJn{t) = —U{t,z) 

z=0 

then, as follows from derivation of F(j, t; cj), the quantity 

F(j,t;u;) = D(j,t)exp|-^^ dt{uj{t, z) - l)"^ (12) 
will describe the pre-confinement partons multiplicity distribution in the frame of LLA. 



Inserting ([121) into the integral ([ToD , one can find that, if we will use the designation: 

a;(t, z) = dt uj{t, z), 



J to 

j — 1 = {iu{t, z) / ln{l/ x)Y^'^ are essential in it. So, to justify the LLA approximation, one 
should assume that the "mobility" 

{\n{l/x)/iu{t,z)}:^l (13) 

decreases with z or, it is the same, with the multiplicity n. This is the reason why the LLA 
for considered DIS kinematics has a restricted range of validity in the VHM region. 

Nevertheless, in the frame of LLA conditions, as follows from (|l2|) , the generating func- 
tional J-'ab{x, t; z) has the following estimationdn ^^^^(x, t; z) oc {ln{l/x)u{t,z)y/^. Therefore, 
since the coupling is a constant, InTabix^i] z = 1) = hi'Dab{x,t) oc (t ln(l/a;))^/^. This is a 
well known result. 

4. The higher correlation functions N^l^'^'"^''{ti,t2, .■.,t^;x,q'^) must be taken into ac- 
count if z is large: the smallness of this correlation functions can be compensated by the 
large value of uj^{t,z). This leads to dominance of the multi-jet processes with jets of ap- 
proximately equal masses. Such kinematics is outside of the LLA. Let us consider it in 
detail. 

It is known H] that u{t, z) is singular at 

Zs-lr^ l/%(t), (14) 

where Ujit) is the mean multiplicity in the jet: Innj(t) ~ \/t. As follows from (|p, the 
singularity ([T^ gives the estimation: = l/nj{t) + 0{l/n), i.e. the character of singularity 
is not important with 0{l/n) accuracy. 

This means that, for instance, the multi-jet event with k > 1 gives the same contribution 
to /i'. The difference may appear only in the pre-exponential factor. But actually the 
difference will appear since the singularity moves with t, Zs = Zs{t). Indeed, if the energy 
conservation law is taken into account and if the produced k jets have the same energy, 
then nf\t) < nj{t). Consequently, zi''){t) > Zs{t) and /j.'j'^^ > n'j. This means that the 
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heavy jet productions mechanism would dominate in the VHM region in accordance with 
the proposition (I). 

5. The above derived result means the increase of the mean transverse energy with 
the multiplicity n. The continuation of this result into the asymptotics over multiplicity 
leads to the deduction that only two heavy jets with masses ~ Vs/2 are produced in the 
deep asymptotics over n. However, the value of such a process would be extremely small, 
oc {as{s)/s). 

Therefore, according to (I) it is most probable to expect the increase of the transverse 
energy with multiplicity, the experimental confirmation of this result was given in [0] , but 
apparently the limiting case of two-jet kinematics will never be exceeded. 
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